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The phase diagram for the vortex states of high-T c superconductors with point defects in B || c 
axis is drawn by large-scale Monte Carlo simulations. The vortex slush (VS) phase is found between 
the vortex glass (VG) and vortex liquid (VL) phases. The first-order transition between this novel 
normal phase and the VL phase is characterized by a sharp jump of the density of dislocations. The 
first-order transition between the Bragg glass (BG) and VG or VS phases is also clarified. These 
two transitions are compared with the melting transition between the BG and VL phases. 
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Vortex states in high-T c superconductors in B \\ c 
axis have been intensively studied. Although the melting 
transition in pure systems has now been understood very 
well, experimental phase diagrams are more complicated 
owing to effects of impurities. In the present Letter, point 
defects are taken into account as the simplest impurity. 

A schematic phase diagram of high-T c vortex states 
with point defects is given in the inset of Fig. [j]. 
These three phases have been observed in YBa2Cu307_5 
(YBCO) @,| and Bi 2 Sr 2 CaCu 2 8+2y (BSCCO). [||] 
The Bragg glass (BG) phase || is characterized by the 
power-law decay of correlation functions of vortex posi- 
tions HH and the triangular Bragg pattern of the struc- 
ture factor. The vortex glass (VG) phase was first defined 
on the basis of phase variables, |?J and can be defined al- 
ternatively on the basis of vortex positions. || Recent 
simulations including screening effects suggest the 

absence of the phase-coherent VG. However, the posi- 
tional VG is expected to be more stable than the phase- 
coherent VG. H Therefore, whether the VG exists as the 
thermodynamic phase or not is still unsettled. 

The BG-VG phase boundary was studied analytically, 
[|iT| -p^| essentially based on the Lindemann criterion. 
Numerically, difference between the BG and VG phases 
was discussed, 0-16| and the overall phase diagram was 
obtained recently. [17jl8| However, studies based on ther- 
modynamic quantities are still lacking. 

Quite recently, another phase has been observed ex- 
perimentally. Nishizaki et al. have found a sharp jump 
of the resistivity |l9| , po[| and the magnetization (2(J in 
optimally-doped YBCO in the vortex liquid (VL) region. 
They pointed out that this might be the transition to the 
vortex slush (VS) phase originally found in irradiated 
YBCO. The VS phase is characterized by evolution of a 
short-range order, and the boundary between the VS and 
VL phases terminates at a critical point. Similar anoma- 
lies were also observed in BSCCO. (2§j2| 

In the present Letter, we compose the phase diagram 
in the pinning-strength (e)-temperature (T) plane (Fig. 
|l|) on the basis of thermodynamic quantities. The first- 
order melting transition occurs between the BG and VL 
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FIG. 1. e-T phase diagram of the model with r = 20. 
Shaded phases exhibit superconductivity. Schematic B—T 
phase diagram is shown in the inset for comparison. 



phases as in pure systems. p6| , |27|l A first-order transition 
line stretches from the melting line into the VL region 
dividing the VS and VL phases. The VG phase exists 
at much lower temperatures. The boundary between the 
BG and VG or VS phases is almost independent of tem- 
perature, and the phase transition is of first order. This 
phase diagram is quite similar to the B—T phase diagram 
observed experimentally. 

Our formulation is based on the three-dimensional 
frustrated XY model on a simple cubic lattice, |28|j2rj|] 
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with the periodic boundary condition. Screening effects 
are not included in this model. A uniform magnetic field 
B is applied along the c axis, and the averaged number 
of flux lines per plaquette is given by / = l 2 B/$>o. Here / 
stands for the grid spacing in the ab plane. Point defects 
are introduced as the plaquettes which consist of four 
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weaker couplings and are randomly distributed in the ab 
plane with probability p. Couplings are given by Jij = 
(1 — e) J (0 < e < 1) on the point defects, and = J 
elsewhere. Here we concentrate on the model with T = 
20, f = 1/25 and p = 0.003. We do not consider the 
lower critical point [p9j30| at present. The system size 



is L x = L y = 50 and L c = 40, which is large enough to 
analyze the pure system (e = 0). 

For each e, Monte Carlo simulations are started from a 
very high temperature, and systems are gradually cooled 
down. After such annealing, further equilibration and 
measurement are performed at each temperature. Typi- 
cal simulation times are 4 ~ 12 x 10 T Monte Carlo steps 
(MCS) for equilibration, and 2 - 5 x 10 7 MCS for mea- 
surement. The present simulations are based on one sam- 
ple. Since configurations of point defects are independent 
in different ab planes, it is reasonable to expect that there 
exists a self- averaging effect. We calculate the internal 
energy e, the specific heat C, the helicity modulus along 
the c axis T c , [|28|]26[ and the phase difference between 
the nearest-neighbor ab planes (cos(</> n — </> n+ i)), together 
with the ratio of entangled flux lines to total flux lines 
^Vent/^Vflux, H?]] the density of dislocations in the ab plane 
Pd, and the structure factor of flux lines in the ab plane. 
The helicity modulus is proportional to the superfluid 
density, and is the order parameter of superconductivity. 
The inter-layer phase difference is related [|3l] to the fre- 
quency of the Josephson plasma resonance (JPR). [ f32p3[ 

BG-VL and VS-VL phase transitions. — Temperature 
dependence of e, C, T c and (cos(</> n — n+ i)) at e = 0.05, 
0.07 and 0.11 is displayed in Figs. |(a)-(d). The BG-VL 
transition is observed at e = 0.05, and the VS-VL one 
at e = 0.07 and 0.11. Jumps of e and (cos(</> n — (/> n +i)) 
and the (^-function peak of C are observed at the melting 
temperature T m or the slush temperature T 8 \. Finite la- 
tent heats Q = Ae indicate that these phase transitions 
are of first order. The anomalies of the VS-VL transition 
at e = 0.07 are as large as those of the BG-VL transition 
at e = 0.05, and gradually lose sharpness as e is increased 
(see Figs. |(c) and |(d)). At e = 0.14, no anomalies are 
observed other than a small jump of (cos(0 n — n+ i)). 
These properties can be explained by the existence of 
the critical point Jl9|— which terminates the first-order 
VS-VL transition line around e = 0.14. 

The most important difference between the BG-VL 
and VS-VL transitions is seen in T c . In the BG-VL 
transition, this quantity appears discontinuously at T m . 
On the other hand, T c remains vanishing for T < T s \ in 
the VS-VL transition. In other words, the BG-VL phase 
transition is the superconducting-normal one, while the 
VS-VL transition occurs between two normal phases. 
The latter transition does not contradict the existence of 
the critical point. The proliferation of T c at lower tem- 
peratures signals the phase transition to the VG phase, 
though error bars of T c are fairly large owing to very 
large correlation time in the VG phase. 
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FIG. 2. Temperature dependence of (a) e at e = 0.05 (x), 
0.07 (+) and 0.11 (o), and C (A), T c (□) and (cos(0 n -0 n+ i)) 
(o) at (b) e = 0.05, (c) 0.07 and (d) 0.11. The origin of e is 
varied for each e in (a). 
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Direct observation of flux lines also reveals the differ- 
ence between these two first-order transitions. Temper- 
ature dependence of N ent /Nfi UK and pa at e = 0.05 and 
e = 0.11 is displayed in Figs. [|(a) and ||(b), respectively. 
^Vent/^Vflux shows a sharp jump at T m (Fig. ||(a)) as in 
pure systems, |27j while its temperature dependence is 
continuous around T s \ (Fig. ||(b)). The quantity pa shows 
sharp jumps both at T m and T 8 \ as predicted by Kierfeld 
and Vinokur. Q These properties are consistent with 
the structure factors shown in the same figures. A ring- 
like pattern is seen in the VL phase as in pure systems 
[ p6[ both at e = 0.05 and 0.11. The clear triangular Bragg 
pattern for T < T m at e = 0.05 represents the formation 
of a hexatic quasi long-range order. The obscure Bragg 
pattern with a 6-fold symmetry for T < T s \ at e = 0.11 
might stand for domains of a short-range hexatic order 
in the ab plane divided by dislocations. 

BG-VG and BG-VS phase transitions. — Pinning- 
strength dependence of e and (cos(0 n — </> n +i)) is shown 
in Figs. |(a) and |(b) for T = 0.06, 0.07 and 0.08J/&B- 
Finite latent heats Q indicate that these phase transitions 
are of first order. The quantity (cos(</> n — n+ i)) jumps 
simultaneously, which is consistent with Gaifullin et a/.'s 
JPR experiment of BSCCO. |35| Since the anomalies are 
always observed between e = 0.065 and 0.070, the phase 
boundary is almost independent of temperature as shown 
in Fig. |l|. We also find a sudden jump of AT ent /iVfl ux on 
this phase boundary as predicted theoretically. [p"l|jL3[ 

The jumps of e and (cos(</> n — n+ i)) (abbreviated as 
A (cos) from now on) on this almost-flat phase boundary 
are about one order smaller than those on the melting 
line (see Figs. ||(b) and ||(b)). The ratio of the jump of 
the Josephson energy Aej = — (J/r 2 )A(cos) to Q reveals 
the difference between these two transitions clearly. As 
in extremely anisotropic pure systems, J3(J this ratio is 
about one half on the melting line, e.g. Aej « 5.5 x 10 -4 J 
and Q « 1.1 x 10" 3 J at e = 0.05 (see Figs. |(a) and 
||(b)). On the other hand, Aej/Q is about unity on this 
almost-flat phase boundary, e.g. Q « Aej « 9.0 x 10 -5 J 
at T = 0.07 J/k B (see Figs. |(a) and |(b)). The latter 
result shows a sharp contrast to the JPR experiment of 
BSCCO, where Aej/Q ^> 1 seems to be satisfied. [ |37| 

Discussions. — On the melting line of pure systems, 
A (cos) is proportional to T 2 f and gradually approaches 
a saturated value, f§ A(cos) w 0.3. [[||35| When the 
anisotropy is as small as that of YBCO (r « 7 ~ 8), 
A (cos) is small both on the melting line and on the BG- 
VG phase boundary. In the present system (r = 20), 
A(cos) (~ 0.22 at e = 0.05) is as large as the sat- 
urated value on the melting line, while it is small on 
the BG-VG/VS phase boundary. Therefore, a jump in 
(cos((/> n — (j) n+ i)) inevitably occurs in the VL region, 
which results in the VS-VL phase transition. On the 
other hand, when the anisotropy is as large as that of 
BSCCO (r > 150), A (cos) has reached the saturated 
value both on the melting line and on the BG-VG phase 
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FIG. 3. Temperature dependence of A^ent/^Vfl ux (o) and 
pa (•) at (a) e = 0.05 and (b) e = 0.11. Structure factors at 
T = 0.087 and 0.088 J/kB are displayed in (a), and those at 
T = 0.080 and 0.083 J/k B are displayed in (b). 
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FIG. 4. Pinning-strength dependence of (a) e and (b) 
(cos(0 n -0 n+ i)) at T = 0.06 (□), 0.07 (o) and 0.08J/fe (A). 
The origin of e is varied for each T in (a). 
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boundary as shown experimentally. [[35| In such a case, it 
might be difficult to observe a jump of (cos(0 n — n+ i)) 
outside of the BG phase by the JPR. [p9| , |37[ 

Finally, the present results summarized in Fig. [l] are 
compared with theoretical studies related to the VS 
phase in literature. When Worthington et al. (2l) pro- 
posed the VS-VL transition line as a reminiscent of the 
melting line in pure systems, the BG phase was out of 
the scope. Ikeda Q derived a phase diagram consist- 
ing of the VG, VS and VL phases, and argued that the 
BG phase and the VS phase cannot coexist. Quite re- 
cently, he modified his argument and proposed a possi- 
ble phase diagram including both the BG and VS phases. 
|[IH However, he simply assumed the existence of the BG 
phase in this article. Kierfeld and Vinokur [|2| obtained 
a phase diagram consisting of the BG phase and a first- 
order transition line stretching from the melting line into 
the VL region with a critical point. Although they inter- 
preted this transition line as the VG-VL one, it turns out 
to correspond to the VS-VL one as shown in the present 
study. Reichhardt et al Q numerically found a window- 
glass-like region with diverging time scales and a finite 
correlation length, which might also correspond to the 
VS phase. We believe that the present simulations are 
the first theoretical derivation of the phase diagram con- 
sisting of four phases (Fig. ||) based on thermodynamic 
quantities and from a single model. Although finite-size 
analyses of T c have not been performed at present, 
our system size, L c = 40 for T = 20 and / = 1/25, is 
already very large. According to the anisotropy scaling, 
[Q this size corresponds to L c « 113 for the parameters 
in Ref. @, T = and / = 1/20. 
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